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Abstract 

Classifications of twin primes are established and then applied to 
triplets that generalize to all higher multiplets. Mersenne and Fer- 
mat twins and triplets are treated in this framework. Regular prime 
number multiplets are related to quadratic and cubic prime number 
generating polynomials. 
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1 Introduction 



There exist extensive tables of twin, triplet and quartet primes. There are 
no such systematic analyses of higher multiplets. Here we outline a more 
systematic analysis of generalized twin primes, triplets and regular multiplets 
that are connected with prime producing quadratic and cubic polynomials. 
We follow standard practice ignoring as trivial the prime pairs (2,p) of odd 
distance p — 2 with p any odd prime. 
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2 Classifications of twin and triplet primes 



Definition 2.1. A triplet Pi } p m = Pi + 2di,p/ = p m + 2g?2 of odd prime 
numbers with pi < p m < pf is called a generalized triplet. 

Each triplet consists of three generalized twin primes (pi,p m ), (p m ,Pf), 
(pi,Pf). Empirical laws governing triplets therefore are intimately tied to 
those of the generalized twin primes. 

2.1 Twin prime classifications 

There are two schemes of parametrizations of twin primes which we now 
describe. 

Theorem 2.2. Let 2D be the distance between odd prime numbers Pi,Pf 
of the pair, D a natural number. Then there are three mutually exlusive 
classes of generalized twin primes that are parametrized as follows. 

I: Pi = 2a- D, p f = 2a + D, D odd; (1) 
II: Pi = 3(2a- 1) - D, p f = 3(2a - 1) + D, 2\D, 3 J(D; (2) 
III: Pi = 2a + l-D, p f = 2a + l + D, 6\D, D > 6, (3) 

where a is the running integer variable. Values of a for which a prime pair 
of distance 2D is reached are unpredictable (called arithmetic chaos). 

Each of these three classes of this classification [1] contains infinitely many 
(possibly empty) subsets of prime pairs at given even distance. 

Proof. Let us first consider the case of odd D. Then p^ = 2a — D for some 
positive integer a and, therefore, pf = Pi + 2D = 2a + D. The median 2a is 
the running integer variable of this class /. 

For even D with D not divisible by 3, let pi = 2n + 1 — D so that 
Pf = 2n + 1 + D for an appropriate integer n. Let pi ^ 3, thus excluding a 
possible first pair with pi = 3 as special. Since of three odd natural numbers 
at distance D from each other one is divisible by 3, the median 2n + 1 must 
be so, hence 2n + l = 3(2a — 1) for an appropriate integer a. Thus, the median 
3 (2a — 1) of the pair 3 (2a — 1) ± D is again a linear function of a running 
integer variable a. These prime number pairs constitute the 2nd class //. 

This argument is not valid for prime number pairs with 6\D, but they 
can obviously be parametrized as 2a + 1 ± 6d, D = 6d. They comprise the 
3rd and last class III of generalized twins. Obviously, these three classes are 
mutually exclusive and complete except for the special cases, o 
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Example 1. 

Ordinary twins 2a ± 1 for a = 2, 3, 6, 9, 15, . . . have D — 1 and are in class 
/. For £> = 3, 2a ± 3 are twins for a = 4, 5, 7, 8, 10, ... . For D = 5, 2a ± 5 
are twins for a = 4, 6, 9, 12, 18, ... . No twins are ever missed or special in 
this class /, an advantage of this classification. 

For D = 2, 3(2a - 1) ± 2 are twins in class II for a = 2, 3, 4, 7, 8, ... . 

For D = 6, 2a + 1 ± 6 are twins in class J/J for a = 5, 6, 8, 11, 12, ... . 

Special twins are 5 ± 2, 7 ± 4, 11 ± 8, . . . 

The second classification of generalized twins involves arithmetic progres- 
sions of conductor 6 as their regular feature. It is well known that, except 
for the first pair 3, 5, ordinary twins all have the form (6m — 1, 6m + 1) for 
some natural number m. They belong to class /. 

Example 2. 

Prime pairs at distance 2D = 4 are in class II and of the form 6m + 
1, 6(m+l) — 1 for m = 1, 2, 3, . . . except for the singlet 3, 7. At distance 2D = 6 
they are in class / and have the form 6m — 1, 6(m + 1) — 1 for m = 1, 2, 3, . . . 
that are intertwined with 6m + 1, 6(m+ 1) + 1 for m — 1, 2, 5, .... At distance 
2D = 12 they are in class /// and of the form 6m — l,6(m + 2) — 1 for 
m = 1, 3, 5, . . . intertwined with 6m + 1, 6(m + 2) + 1 for m = 1, 3, 5, . . . . 

In general, the rules governing the form 6m ± 1, 6m + b depend on the 
arithmetic of D and a making the second classification of generalized twin 
primes difficult to deal with generally. We now apply Theor. 2.2. 

Theorem 2.3. Let 2D be the distance between odd prime numbers Pi,Pf 
of the pair. Then for class III, D = (mod 6) and pf^ = 2a + 1 (mod 6). 
If a = (mod 3) then pf^ = 1 (mod 6). If a = — 1 (mod 3) then pf^ 
I (mod 6). 

For class II, D = 2r (mod 6), r = ±1 yields Pi = ±1 (mod 6),Pf 
'■ I (mod 6); r = ±2 gives pi = =pl, pf = ±1 (mod 6). 

For class I, and D = 1 + 2r (mod 6), r = 0, ±1, a = a (mod 3) the 
prime pair obeys pi = — 1 — 2r — 2a (mod 6),p/ = 1 + 2r + 2a (mod 6). 
For r = and r = —1, i.e. D = ±1 (mod 6), ao = 0. For r = 1, ie. 
.D = 3 (mod 6), a ^ (mod 3). 

Proof. For class a = 1 (mod 3) is ruled out because then pfi = 3 
(mod 6). For class //, = 3(2a — 1) ± D, so D = 2r (mod 6) implies 
j = — 3 ± 2r (mod 6) for r = ±1, etc. For class /, and D = 1 (mod 6) 
the cases a = ±1 (mod 3) are ruled out because they imply either 3\pi or 
3\pf, except a = 2, i.e. D = 2. This also is the case for D = — 1 (mod 6). 
For D = 3 (mod 6) the case a = (mod 3) is obviously ruled out. o 
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This concludes the classifications of twin primes at constant distances. 
2.2 Triplet prime classifications 

Generalized triplets have the form 6m ±1, 6m + ai, 6m + a 2 , except for singlet 
exceptions and this generalizes to prime quadruplets, quintuplets, etc. 

Rules for singlets or exceptions among generalized triplet primes are the 
following. 

Theorem 2.4. (i) There is at most one generalized prime number triplet 
with distances [2D, 2D] for D = 1, 2, 4, 5, . . . and 3 J(D. 

(ii) When the distances are \2d\,2d 2 \ with 3\d 2 — d\, and 3 j{d\, there is 
at most one triplet p { = 3,p m = 3 + 2d±,pf = 3 + 2c?! + 2d 2 for appropriate 
integers d\,d 2 . 

Proof, (i) Because of three odd numbers in a row one is divisible by 
3, 3, 5, 7 is the only triplet at distances [2, 2] and, for the same reason, there 
is only one triplet 3, 7, 11 at distances [4,4], one only at distances [8,8] i.e. 
3, 11, 19, at [10, 10] i.e. 3, 13, 23 and in general at distances [2D, 2D] for D 
not divisible by 3. The argument fails when 3\D. (ii) Of Pi, p m = Pi + 2di 
(mod 3), pf = pi + Adi (mod 3) at least one is divisible by 3, which must 
be pi. o 

Naturally, the question arises: Are there infinitely many such singlets, i.e. 
exceptional triplets? 
Example 3. 

At distances [2, 8], the singlet is 3, 5, 13 and at [8, 2], it is 3, 11, 13. 

At distances [2,4], the triplets are 2n — 3,2n — l,2n + 3 with n = 1 
(mod 3) and at [4, 2], they are 2n — 3, 2n+ 1, 2n + 3 with n = — 1 (mod 3). 
Writing n = 3m ± 1 for these cases, we obtain 6m — 1, 6m + 1, 6m + 5 and 
6m — 5,6m — 1,6m + 1 for these triplets, respectively. At distances [2,4], 
triplets occur for n — 4, 7, 10, . . . i.e. m — 1,2,3, ... , while at [4, 2], they are 

at n — 5, 8, 20, . . . i.e. m — 2,3,7, These triplets are in the classes (/, II) 

and (II, I), respectively, with I, II denoting symbolically the classes of the 
first generalized twin prime classification. 

At distances [2, 6] the generalized triplets are 2n — 3, 2n — 1, 2n + 5 with 
n = 1 (mod 3) or 6m — 1, 6m + 1, 6m + 7. The triplet 3, 5, 11 for n — 3 is 
the only exception. For [6, 2] they are 2n — 5, 2n + 1, 2n + 3 with n = — 1 
(mod 3) or 6m — 1, 6m + 5, 6m + 7. They all are in the class (J, /). 

Applying Theors. 2.2, 2.3 we obtain the following triplet classifications. 



4 



Corollary 2.5. (i) The class (1,1) is made up of the triplets 

Pi = 2a- D ± , p m = 2a + D 1 = 2b- D 2 , Pf = 2b + D 2 (4) 

with odd Di, D 2 and a, b appropriate integers subject to b — a = (D\ + D 2 )/2. 
Hence 

p f = 2a + D 1 + 2D 2 , (5) 

and the prime number pair 

( Pf , Pi )=2a + D 2 ±(D 1 + D 2 ) (6) 

is in class II, or special, or III. 
(n) U 

£>i = 2ri + 1 (mod 6), n = 0, ±1; D 2 = 2r 2 + 1 (mod 6), r 2 = 0, ±1, 
a = a (mod 3), a = 0, ±1 (7) 

Pi = — 2r\ — 1 + 2a (mod 6), p m = 2r\ + 1 + 2a (mod 6), 

p f = 2a + 3 + 2ri + 2r 2 (mod 6), D 1 + D 2 = 2(r\ + r 2 ) + 2 (mod 6). 

(8) 

Now specific cases (ii) can be worked out by substituting values for r\, a : 
D 1 = 1 (mod 6), a = (mod 3) yield Pi = —1 (mod 6), p m = 1 
(mod 6), pf = 1 + 2D 2 (mod 6) = ±1 (mod 6) with D 2 = 1 (mod 6) 
ruled out, etc. Of course, Example 3 is consistent with this. 

Corollary 2.6. (i) The class (I, II) consists of prime triplets at distances 
Di odd and D 2 even so that 

Pl = 2a- D x , p m = 2a + D 1 = 3(26 - 1) - D 2 , p f = 3(26 - 1) + D 2 (9) 

with appropriate a, b subject to a = 3b — \(D X + D 2 + 3). Hence pf = 2a + 
Di + 2D 2 and the twin 

( Pf , Pi ) = 2a + D 2 ±(D 1 + D 2 ) (10) 

is in class I. 
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(n) U 



D 1 = 2n + 1 (mod 6), n = 0, ±1; D 2 = 2r 2 (mod 6), r 2 = ±1, 
a = a (mod 3), a = 0, ±1 (11) 



then 



Pi = —2r\ — 1 + 2a (mod 6), p m = 2r\ + 1 + 2a (mod 6) 

= -3-2r 2 (mod 6), p f = -3 + 2r 2 (mod 6). (12) 

Corollary 2.7. Tne c/ass (/J, //) has Di even and D 2 even with both 
Di not divisible by 3. The triplets are 

Pl = 3(2a - 1) - D u Pm = 3(2a - 1) + D 1 = 3(26 - 1) - D 2 , 

Pf = 3(2b-l) + D 2 (13) 

so that3(b-a) = |(Di + D 2 ). Therefore 3|Di + D 2 and Di = -D 2 (mod 3). 

P/ = ^ + 2(Di + D 2 ) = 3(2a - 1) + D l + 2D 2 , 
(p f , Pi ) = 3(2a-l) + D 2 ±(D 1 + D 2 ) (14) 

is in class III with 6\Di + D 2 . 
(n) U 

Di = 2r-i (mod 6), n = ±1 (15) 

then 

Pi = — 3 — 2ri (mod 6), p m = — 3 + 2r x (mod 6) = — 3 — 2r 2 (mod 6), 
p/ = — 3 + 2r 2 (mod 6) = ^ (mod 6). (16) 

Corollary 2.8. The class (I, III) has the triplets 

p t = 2a- D 1 , p m = 2a + D l = 2b + 1 - 6d 2 , p f = 2b + 1 + 6d 2 (17) 

where 2b + 1 = 2a + D x + 6d 2 . So p/ = 2a + D x + 12d 2 and 

(p f ,p t ) = 2a + Qd 2 ±(D 1 + 6d 2 ) (18) 

is m c/ass /. 
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(H) U 

D 1 = l + 2r 1 (mod 6), ri = 0,±l;a = ao (mod 3) (19) 

then 

Pi = —1 — 2r\ + 2a (mod 6), p m = 26 + 1 (mod 6) 

I + 2r\ + 2a (mod 6), pf = p m (mod 6), 2b = 2{r\ + a ) (mod 6). 

(20) 

Corollary 2.9. C/ass (77, 777) has D\ even and not divisible by 3 and 
D 2 = 6d 2 with the triplet form 

Pi = 3(2a -!)-£>!, p ro = 3(2a - 1) + £>i = 26 + 1 - 6d 2 , 

p f = 26 + 1 + 64 (21) 

so that 

2b + l = 3(2a - 1) + D 1 + 6d 2 , p/ = 3(2a - 1) + D x + 12d 2 . (22) 

So 

(p f , p^ = 3(2a - 1) + Qd 2 ± (Dx + Qd 2 ) (23) 

«s m class II. 
(H) If 

D 1 = 2r 1 (mod Q),r 1 = ±1, a = a (mod 3), a = 0, ±1 (24) 

then 

Pi = — 3 — 2r\ (mod 6), p m = — 3 + 2r\ (mod 6) 

= 26+1 (mod 6), p f = p m (mod 6). (25) 

Corollary 2.10. The class (777,777) has the triplet form 

Pi = 2a + 1 - 6di, p m = 2a + 1 + 6di = 26 + 1 - 6d 2 , 26 + 1 + 6d 2 (26) 

/or appropriate a, 6 so £/ia£ 

6 = a + 3(di + d 2 ), P/ = 2a + 1 + 6(d 1 + 2d 2 ) (27) 
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and 



(PfiPi) — 2a + 1 + 6d 2 ± 6(di + d 2 ) 



(28) 



is in class III, too. 

(ii) If a = a (mod 3), b = b (mod 3) then 



Pi = 2a + 1 (mod 6) = p m (mod 6) = pf (mod 6). 



(29) 



The classes (//, J), (///, I), (III, II) are handled similarly. Several examples 
have been given above. These nine classes of generalized prime number 
triplets are mutually exclusive and complete except for the singlets. 

These twin and triplet prime classifications represent regularities that 
generalize to quadruplet primes which come in 3 3 mutually exlusive classes, 
quintuplet primes in 3 4 such classes except for singlets, etc. 

3 Special twin and triplet primes 

Here we consider Mersenne and Fermat twins and triplets. 
3.1 Mersenne twins 

A simple application of the second classification is the following 

Corollary 3.1. If 2 P — 1, with an odd prime number p, is a Mersenne 

prime, then TP + 1 is composite. 

Proof. Since TP — 1 ^ 6m — 1, the pair TP ± 1 is not a twin prime, o 
Of course, it is well known that 3|2 P + 1 but this requires an algebraic 

identity for the factorization: 



Let us now consider Mersenne twins with TP — 1 as the first prime number 
of the pair. 
Example 4- 

The pair 2 P — 1, 2 p +5 is a Mersenne pair in class / for p — 3, 5, . . . , p = — 1 
(mod 3) and p > 7. The qualification is due to the factorization in Lemma 3.2. 

The pair TP — 1, 2 P + 9 is a Mersenne pair in class / for p — 2, 3, 5, 7, 

The pair TP — 1, TP + 3 is a Mersenne pair in class // for p — 2, 3, 7, . . . , 
p = — 1 (mod 4). The restriction is due to Lemma 3.3. 



a p + 1 = ( a + l)(a : 



a p - 2 ±... + a 2 -a + l). 



(30) 
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The pair 2 P — 1, 2 P + 11 is a Mersenne pair in class III for p — 3, 5, 7, . . . . 
Lemma 3.2. 7/p = 1 (mod 3) and p > 7 then 7\2 P + 5. 
Proof. This follows from the factorization 

2 P + 2 3 - 2 - 1 = (2 3 - 1)(2 P ~ 3 + 2 P " 6 + • • • + 2 + 1). o (31) 

Lemma 3.3. If p = 1 (mod 4) then 5|2 P + 3. 
Proof. 

/(p-i)/2 \ 
2 p + 2 2 -2 + l = (2 2 + l) f £ (-l) j_1 2 p - 2 ^ + ll . ❖ (32) 

Next we list Mersenne twins with 2 P — 1 as the second prime number of 
the pair. 

Proposition 3.4. (i) The pair2 p — 5, 2 P — 1 is a Mersenne pair in class II 
forp = 3 only. There are no Mersenne prime twins 2 P — 2 2n+l — 3, 2 P — 1; n — 
1,2,... except p = 3, n = 0. (it) There are no Mersenne twins 2 P — 3, 2 P — 1 
ui/jen p = —1 (mod 4) except for the pair 5, 7. fra) There are no Mersenne 
twins 2 P — 7, 2 P — 1 w/ien p = 1 (mod 4). 

Proof, (i) This holds because 3|2 P — 5 for p > 3 which follows from the 
first factorization 

/p-i 

2 P - 2 - 2 - 1 = (2 + 1) ]T (-2) p - m - 1 

\m=l 

(p-2n-l \ 
(-2) p - m - ll , 

p> 2n + 3, n = 1,2,.... (33) 

and the next cases from the second factorization. The case n = 1 gives 
3|2 p - 11, p > 5. 

(ii) is due to the factorization, for p = —1 (mod 4), 

/(p-1)/2 \ 

2 P - 2 2 + 2 - 1 = (2 2 + 1) I E (-l)™-^^ 2 ™ - 1 j . (34) 

(iii) follows from the factorization, for p = 1 (mod 4), 

/(p-i)/2 \ 
2 f-2 2 -2-l = (2 2 + l) I E (-l) m " 1 2 p - 2m -lj . o (35) 



Example 5. 

The pair 2 P - 15, 2 P - 
The pair 2 P - 19, T> - 
The pair 2 P - 13, 2 P - 
The pair 2? - 25, 2? - 



1 is in class /, and p = 
1 is in class /, and p = 
1 is in class III, and p 
1 is in class III, and p 



5, 7 are such cases. 
5, 7 are such cases. 
= 5, 13 are such cases. 
= 5, 7, 13 are such cases. 



3.2 Fermat twins 

Here we consider Fermat prime pairs with the Fermat prime being its first 
member. 

Example 6. 

2 2 ™ + 1, 2 2 " + 3 are twins in class I for n = 0, 1, 2, 4, ... , i.e. 2 n ~ l ^ 1 
(mod 3), n > 1. The qualification is due to (i) in Lemma 3.6. 

The pair 2 2 ™ + 1, 2 2 ™ + 7 is in class / and n = 1, 2, 3 are such cases. 

The pair 2 2 " + 1, 2 2 " + 13 is in class 77 and n = 1, 2, 3 are such cases. 

Corollary 3.5. (i) The twin prime?), 7, is the only one of'2 2n + 1, 2 2 " + 5 
in class II forn = 0. (ii) The pair 2 2 " + 1, 2 2 " + 9 is in class III and n — 0, 1 
are the only such cases. 

Proof, (i) follows from the factorization 

(2 n -l \ 
E (-l) J_1 2 2n_j + lj , n > 0. (36) 

and (ii) from 

2 2 " + 2 2 + 2 2 + 1 = (2 2 + 1) ^ jjVl)'--^ 2 **- 2 '' + 1 j , n > 0. o (37) 

Lemma 3.6. (%) If 2 n = 2 (mod 3), n > 1 t/ien 7|2 2 " + 3. ^ If 2 n = \ 
(mod 3) toen 7|2 2 " + 5. 

Proof, (i) follows from the factorization 

2 2" + 2 3 - 2 2 - 1 = (2 3 - 1) (2 2 "- 3 + 2 2 "" 6 + . . . + 2 2 + l) , 

2 n = 2 (mod 3), (38) 

and (ii) from 

2 2 " + 2 3 - 2 - 1 = (2 3 - 1) (2 2n ~ 3 + 2 2 "- 6 + . . . + 2 + l) , 

2 n = 1 (mod 3). o (39) 
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Fermat prime pairs with the Fermat prime being its second member are the 
following. 

Example 7. 

The pair 2 2 ™ — 5, 2 2 ™ + 1 is in class I and n = 2, 3 are such cases. 

The pair 2 2 ™ — 3, 2 2 ™ + 1 is in class II and n = 2 is such a case. 

The pair 2 2 ™ — 11, 2 2 ™ + 1 is in class III and n = 2 is such a case. 

Proposition 3.7. There is no Fermat twin primes of the form 2 2 ™ — 
1, 2 2 " + 1; 2 2 " - 7, 2 2 " + 1; 2 2 " - 19, 2 2 " + 1; 2 2 " + 1; 2 2 " - 21. 

Proof. This holds because 5|2 2 " - 1, 3|2 2 " - 7, 2 2 " - 19, 5|2 2 " - 21 which 
is based on the following factorizations: 

2 2" _ i = (2 2 + 1) ^ p\-iy-i 2 2n - 2 i - 1 j , n > 

(2™- 1 -m 

1 < m < 2 n -\ n > 1, 

j2 i-iy-^-i - 1 

m—1,2, o 




3.3 Mersenne triplets 

Here we list prime triplets where the Mersenne prime comes first. 
Example 8. 

2 p - i, 2 p + 3, 2 p + 9 yields triplets for p = 2, 3 : 3, 7, 13; 7, 11, 13. 

2 p - 1, 2 p + 5, 2 p + 9 yields triplets for p = 3, 5 : 7, 13, 17; 31, 37, 41 
and2 p -l, 2 p + 5, 2 P + 11 yields triplets for p = 3, 5:7, 13, 19; 31, 37, 43. 

2P-1, 2 p +3, 2 p +ll yields triplets forp = 3, 7:7, 11, 19; 127, 131, 139. 

Corollary 3.8. (i) 2 P - 1, 2 P + 3, 2 P + 7 for p = 2 yields the only such 
triplet 3, 7, 11. (ii) 2 P — 1, 2 P + 5, 2 P + 7 yields no Mersenne triplets, and 
(Hi) 2 P — 1, 2 P + 3, 2 P + 5 which is in class II, I yields no triplets except for 
p = 3, namely 7, 11, 13. 

Proof, (i) and (ii) follow from the factorization 

2 p + 2 2 + 2 + 1 = (2 + 1) ( E(-2) P ~ J + 1 J , V > 2- (41) 
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(iii) follows from Lemma 3.2 and the factorization 

2 p + 2 3 - 2 2 - 1 = (2 3 - l)(2 p - 3 + 2 p - 6 + --- + 2 2 + l), 

p = -1 (mod 3). o (42) 

Prime triplets where the Mersenne prime comes last are the following. 
Example 9. 

For p = 3 the triplet 2 P — 5, 2 P — 3, 2 P — 1 becomes 3, 5, 7 which is a 
singlet case. 

2 p _ 7 ; 2 p - 3, 2 P - 1 yields 23, 29, 31 for p = 5. 

2 p - 13 ; 2 P - 3, 2 p - 1 yields a triplet for p = 5 : 19, 29, 31. 

2P-19, 2 P -15, 2 P -1 yields triplets for p= 5, 7: 13, 17, 31; 109, 113, 127. 

Corollary 3.9. 2 P — 7, 2 P — 5, 2 P — 1 yields no Mersenne triplets. 2 P — 
1, 2 P + 3, 2 P + 5 yields no triplets except for p = 2. 

Proof. This follows from the factorizations (31), (32), (33), (135]) . o 

Corollary 3.10. The triplet 2 P - 3, 2 P - 1, 2 P + 3 is the only one. 

Proof. This follows from the factorization of 5|2 P + 3 in Eq. (3.3) of 
Lemma 3.3 for p = 1 (mod 4) and 5|2 P — 3 for p = — 1 (mod 4) in 
Prop. 3.4. o 

Finally, prime triplets where the Mersenne prime is in the middle are 
composed by the twins where the Mersenne prime is last followed by twins 
where it comes first. 

3.4 Fermat triplets 

We give a few triplets where the Fermat prime comes first. 
Example 10. 

2 2 " + 1, 2 2 " + 3, 2 2 " + 9 yields 3, 5, 11; 5, 7, 13 for n = 0, 1. 

2 2« + 1; 2 2« + 5^ 2 2" + u yieldg 3; 7) 13 for n = Q 

Corollary 3.11 (%) 2 2 " + 1, 2 2 " +3, 2 2 " + 5/orn = 0, yields 3, 5, 7 which 
is the only case, (ii) 2 2 ™ + 1, 2 2 ™ + 5, 2 2 " + 9 for n = yields 3, 7, 11 as t/ie 
only case. 

Proof, (i) follows from the factorization (i) and (ii) from (ii) in Prop 3.7. o 
We list a few triplets where the Fermat prime comes last. 
Example 11. 

2 2 " - 5, 2 2 " - 3, 2 2 " + 1 yields 11, 13, 17 for n = 2. 
2 2 " - 9, 2 2 " - 3, 2 2 " + 1 yields 7, 13, 17 for n = 2. 
2 2n - 11, 2 2 " - 3, 2 2 " + 1 yields 5, 13, 17 for n = 2. 
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2 2 " - 13, 2 2 " - 3, 2 2 " + 1 yields 5, 13, 17 for n = 2. 
Triplets where the Fermat prime is in the middle are composed by twins 
where the Fermat prime is last followed by twins where it comes first. 

4 Regular prime multiplets 

Extensions of the twin and triplet primes at constant distances to quadru- 
plets, quintets, etc. exist but are too numerous to be analyzed systematically 
here. 

4.1 Regular multiplets from quadratic polynomials 

We restrict our attention to those with regularly increasing (or decreasing) 
distances, such as 2N, N = 1, 2, . . . i.e. pi,p2 = Pi+2,p^ = Pi+6, ■ ■ ■ ,Pn+i = 
Pi+n(n+l), . . . , N with the sequence of differences Apj = Pj+i—pj = 2j, j = 
1, 2, . . . , N. Even regular prime triplets and quadruplets are too numerous to 
be listed. We therefore start with quintuplets in Theor. 4.1 below. 

Example 12. There are at least 14 sextets 11, 13, 17,23,31,41; 17, 19, 
23,29,37,47; 41,43,47,53,61,71; 1277,1279,1283,1289,1297,1307; 1607, 
1609, 1613, 1619, 1627, 1637; 3527, 3529, 3533, 3539, 3547, 3557; 21557, 21559, 
21563, 21569, 21577, 21587; 28277, 28279, 28283, 28289, 28297, 28307; 
31247, 31249, 31253, 31259, 31267, 31277; 33617, 33619, 33623, 33629, 33637, 
33647; 55661, . . . , 55691; 113147, 113149, 113153, 113159, 113167, 113177; 
128981, . . . , 129011; 548831, . . . , 548861; 566537, . . . , 566567; 
seven septets 11, ... , 41, 53; 17, . . . , 59; 41, ... , 83; 1277, . . . , 1319; 21577, 
. . . , 21599; 28277, . . . , 28319; 113147, . . . , 113189; 128981, . . . , 129023; five 
octets 11, ... , 53, 67; 17, ... , 73; 41, ... , 97; 21557, . . . , 21599, 21613; 
128981, . . . , 129037; three nonets 11, ... , 67, 83; 17, ... , 89; 41, . . . , 113 
and three decuplets 11, . . . , 83, 101; 17, ... , 117; 41, ... , 131. 
For N = 15 the sequence 17, 19, 23, 29, 37, 47, 59, 73, 89, 107, 127, 149, 173, 
199, 227, 257 is a regular 16-plet. The second such 16-plet starts soon after 
with 41, 43, 47, 53, ... , 281 = 41 + 15-16. In fact, this one extends much longer 
and is the first 40-plet ending with 1523, 1601 = 41 + 39 • 40. There are also 
many almost-regular prime multiplets where just one member is missing, e.g. 
n = 0. 

Naturally, questions arise: Are there infinitely many of these long regular 
multiplets or even longer ones? Are there other long prime multiplets with- 
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out the regular structure imposed by a (quadratic) polynomial on its first 
multiplet which then continues through all of them? Needless to say, the 
existence of so many regular prime multiplets linking primes with each other 
in interlocking ways belies the probabilistic independence of prime numbers 
underlying many conjectures [3], [I] about them. It may not be unreasonable 
to expect that most multiplets repeat infinitely often. Since they are in- 
terlocked their asymptotic distributions are not independent. This suggests 
that asymptotic laws for prime multiplets differ fundamentally from ordinary 
prime numbers. 

The long regular multiplets of Example 12 are related to Euler's prime 
number generating polynomials which, in turn, are related to imaginary 
quadratic number fields over the rationals. There are corresponding polyno- 
mials whose values form regular multiplets that are related to real quadratic 
number fields. Although some of these polynomials have been known for a 
long time with their large number of prime values as the main focus, their 
regular distances within the coherent structure of a regular prime number 
multiplet seem not to have been noted (in print). Except for the multiplet 
aspects many details below are known and documented in Ref. [2J, but some 
are new. 

Theorem 4.1. (i) The Euler polynomials E p (x) = x 2 + x + p with the 
prime numbers p = 2, 3, 5, 11, 17, 41 assume prime number values p + x(x+ 1) 
at distances 2(x + 1) for x = 0, 1, . . .p — 2 forming a regular (p — l)—plet. 

(ii) The polynomials f p (x) = 2x 2 +p with the primes p = 3,5,11,29 
assume prime number values at x — 0, 1, . . . ,p — 1 forming a regular p—plet 
at distances 2{2x + 1). 

(Hi) The polynomials F p (x) = 2x 2 + 2x + |(p + 1) with primes p = 
5, 13,37 = 1 (mod 4) assume prime values for x = 0,1, (p — 3)/2 which 
form a regular (p — i)/2—plet at distances A(x + 1). 

(iv) The polynomials G p , with primes p, q, 

G p (x) = px 2 + px + -{p + q), p<q,pq = 3 (mod 4), 

(p, q) = (3, 5), (3, 17), (3, 41), (3, 89), (5, 7), (5, 23), (5, 47), (7, 13), (7, 61), 
(11, 17), (13, 31) (43) 

assume prime values for x = 0, 1, . . . \{p+q)— 2. which form regular multiplets 
at distances 2p(x + 1) that are independent of the prime q. 

(v) Polynomials corresponding to real quadratic number fields are ga{x) = 
—x 2 + x + |(d — 1) with d > 0, d = 1 (mod 4) and square-free. Their values 
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are prime numbers for x = 2, 3, . . . < \yd — 1 forming regular multiplets at 
distances 2p(x + 1). Relevant values are d = 37, 53, 77, 101, 173, 197, 293, 437, 
677. 

(vi) Quintets, sixtets, septets, octets, nonets and decuplets are generated 
by the following polynomials 



Q p (x) 


= x 2 + x + p, x = 0, . . 


. , 4; p = 11, 17, 41, 347, 641, 1427, 4001, 




4637, 4931, 19421, 22271, 23471, 26711, 27941, 28277, 




31247,32057,33617; 


113147...; 


SX p (x) 


= x 2 + x + p, x = 0, . . 


. 5; p = 11, 17, 41, 1277, 1607, 3527, 




28277, 31247, 33617, 55661, 113147, 128981, 548831, 566537, . . . ; 


S p (x) 


= x 2 + x + p, x = 0, . . 


.,6; p= 11,17,41,1277,28277, 




113147,128981,...; 




O p (x) 


= x 2 + x + p, x = 0, . . 


.,7; p= 11,17,41,128981,...; 


N p (x) 


= x 2 + x + p, x = 0, . . 


.,8; p=ll, 17,41,...; 


D p (x) 


= x 2 + x + p, x = 0, . . 


.,9; p= 11,17,41,.... (44) 



Proof, (i) It is long known [2] that for the listed primes p the values 
Ep(x) = x 2 + x + p at x — 0, 1, 2, ... ,p — 2 are prime numbers. Since 

E p (x + 1) - E p (x) = (x + l) 2 - x 2 + 1 = 2(x + 1), (45) 



the primes p + x(x + 1) form a regular (p — 1)— plet at distances 2(x + 1) 
from each other for x — 0, 1, 2, ... ,p — 2. 

(ii) Since f p (x+l)—f p (x) = 2(2x+l), the primes p+2x 2 , x — 0, 1, . . . ,p— 1 
form a regular p— plet at distances 2(2x + 1). 

(hi) Since F p (x + 1) - F p (x) = 4(x + 1) the primes 2x(x + 1) + (p+ l)/2 
form a regular (p — l)/2— plet at distances 4(x + 1). 

(iv) Since G p (x + 1) — G p (x) = 2p{x + 1) the multiplet structure is clear. 

(v) For g d , this follows from g d (x + l) —gd(x) = —2x. These are multiplets 
with linearly decreasing distances. 

(vi) These multiplets may be verified by a table of prime numbers or 
symbolic-math software, o 

There are many more recently found polynomials in the literature [2] 
which also form regular prime multiplets. 
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4.2 Optimal quadratic polynomials 

This subject has a long history [5] with a rather unsystematic record. It is 
well known that, if P(x) = Sj=o a nX n is a non-constant polynomial of degree 
n > 1 with integral coefficients Oj, [<xq| = a prime number, then P(x) can 
assume prime values at most for x — 0, 1, . . . ,po — 1 because Pol-P(Po)- 

Definition 4.2. The polynomial P(x) = J2 1 j=o a j x '' * s called optimal if 
|P(j)| = pj is prime for j = 0, 1, . . . ,po — 1, forming a regular p — plet. 

Prime number values have to be successive, but they may repeat and be 
negative. This is often caused by negative coefficients in a polynomial. 

As E p (p — 1) = p 2 , the Euler polynomials in (i) of Theor. 4.1 for p = 
2,3,5,11,17,41 are one prime value short of being optimal. But Legen- 
dre's quadratic polynomials [5] f p (x) = 2x 2 + p for p = 3, 5, 11, 29 in (ii) of 
Theor. 4.1 are optimal. 

Since E p (—x) = E p (x — 1), Euler polynomials give repeating prime mul- 
tiplets when they are considered for positive and negative argument. For 
more general polynomials this is not the case. Upon shifting the argument, 
the modified Euler polynomials E p (x — 1) = x(x — 1) +p do become optimal, 
but they repeat the initial prime value. More generally, the identity 

E q {x-n) = x 2 + (1 - 2n)x + E q {n - 1) (46) 

for prime numbers q, p = E q (n — 1) and nonnegative integer n leads to many 
new repeating polynomials with more numerous prime values. For q = 41 
and n = 2, 3, . . . , 40 none of these polynomials is optimal (cf. Eq. (IBTi) ) in 
Cor. 4.5 e.g.), though, including Escott's for n = 40 [2], [5]. 

Definition 4.3. The polynomial P(x) is bi-optimal if |-P(j)| are prime 
for 1 — po < j < po — 1, forming a (2po — 1)— plet. 

This corresponds to combining P(j) for j = 0, 1, . . . ,po — 1 and P{—j) 
for j = 1, . . . ,pq — 1 into one prime number multiplet. For polynomials of 
odd degree the starting value x = is appropriate. For polynomials of even 
degree it is perhaps more natural to include negative arguments as well. 

We now address a question raised by Legendre's f p (x) and Euler 's mod- 
ified polynomials: Are there optimal quadratic polynomials for the missing- 
link primes p = 7, 13, 19, 23, 31, 37? 

For p = 7, Q2(x) = 2x(x — 1) + 7 is optimal, if repeating, generating a 
regular prime number septet with distances AQ 2 (x) = Ax, and it is almost 
bi-optimal because Q2{— x) = Q2^x + 1). 
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For p = 19, Qi{x) = 2x(x — 1) + 19 |5J is not only optimal but almost 
bi-optimal in view of Qj(— x) = Q 7 (x + 1), generating a regular 36-plet with 
the same distance law as Qi- 

For p = 23, Qs(x) = 3x(x — 1) + 23 [5J is not only optimal but almost 
bi-optimal forming a 42-plet because Qs(— x) = Qs(x + 1). 

Proposition 4.4. (i) There is an infinity of optimal quadratic polyno- 
mials for p = 2 : 

Q 1 {x) = x 2 + { Pl -3)x + 2, (47) 

where p\ is prime. If \pi — & \ is prime then Q\(x) is bi-optimal. 
(ii) There are at least three quadratic polynomials forp = 13 : 

Q 3 (x) = x 2 + Tlx + 13, Q 4 (x) = x 2 - 3x + 13, 

Q 5 (x) = 2x 2 -4x + 13, (48) 

that form 12-plets, one prime value short of optimal. Hence 

Q 3 (x - 1) = x 2 + 25x - 13 (49) 

is optimal. Qz{—x) forms a decuplet making Q$(x) almost bi-optimal. There 
is another polynomial, 

Q e (x) = 2x 2 + 26x + 13, (50) 

that forms a decuplet. 

Proof, (i) follows from Qi(0) = 2, Q x (l) = p x and Qi(-l) = 6 - p x . 
(ii) For monic optimal polynomials let 

Q(x) = x 2 + bx+ po, Q(l) =p 1=Po + b + l, Q(2) = p 2 = p + 4 + 2b, (51) 

where Po,Pi,P2 are prime numbers. Then 

b = pi - po - 1, Q = x 2 + (pi - po - l)x + po, p 2 = 2pi - po + 2. (52) 
If Q(x) = 2x 2 + bx + po, then 

Q(l) = px = po + b + 2, Q{2) =p 2 =p + 2b + 8. (53) 

(ii) For po = 13 and p\ = 41 we obtain Q3 from Eq. (1521) and the 12-plet 

Q 3 (l) = 41, Q 3 (2) = 71, Q 3 (3) = 103, Q 3 (4) = 137, Q 3 (5) = 173, 
Q 3 (6) = 211, Q 3 (7) = 251, Q 3 (8) = 293, Q 3 (9) = 337, Q 3 (10) = 383, 
Q 3 (ll) = 431, (54) 
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one prime value short of optimal. Q3 is ascending and non-repeating. 

For p = 13, pi = 11 we get the monic Q 4 , again with a 12-plet of prime 
values 

Q 4 (l) = 11, Q 4 (2) = 11, Q 4 (3) = 13, Q 4 (4) = 17, Q 4 (5) = 23, 
Q 4 (6) = 31, Q 4 (7) = 41, Q 4 (8) = 53, Q 4 (9) = 67, Q 4 (10) = 83, 
Q 4 (ll) = 101, (55) 

which is one prime value short of optimal. It is repeating and non-ascending 
because Q 4 (l) = 11 = Q 4 (2), Q 4 (3) = 13 = Q 4 (0). 

For pi — 11 we get Q§ from Eq. (153]) and the 12-plet 

Q 5 (l) = 11, Q 6 (2) = 13, Q 5 (3) = 19, Q 5 (4) = 29, Q 5 (5) = 43, 
Q 5 (6) = 61, Q 5 (7) = 83, Q 5 (8) = 109, Q 5 (9) = 139, Q 5 (10) = 173, 
Q 5 (ll) = 211, (56) 

which is one prime value short of optimal. It is obviously repeating and 
non-ascending. Since 

Q 5 (-x) = 2x 2 + 4x + 13 = Q 5 (x + 2), (57) 

it generates a decuplet. Therefore, Qs(x) is almost bi-optimal. 
For pi = 41 we get Qq from Eq. (|53|) and the decuplet 

Q 6 (2) = 73, Q 6 (3) = 109, Q 6 (4) = 149, Q 6 (5) = 193, Q 6 (6) = 241, 
Q 6 (7) = 293, Q 6 (8) = 349, Q 6 (9) = 409. (58) 

Corollary 4.5. 27je polynomial 

Q 9 (x) = x 2 + 3x + 19 (59) 

forms a 15-plet for x = 0, . . . , 14, 

Qio(x) = x 2 - x + 11, Qn(x) = 2x 2 + 22x - 11 (60) 

are optimal polynomials for 11, and 

Qu(x) = 2x 2 -4x + 31, Qn{x) = x 2 - 3x + 43 (61) 

/orm 50- and ^2-plets, respectively, one prime value short of optimal. 
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Proof. This follows from Prop. 4.4 in conjunction with 

Q 9 (x-1) = x 2 + x + 17 = E 17 {x), Q 10 (x) = E n {x-l), 
Qii(x) = Q 6 (x - 1), Q 12 (x) = f 29 (x - 1), 
g 13 (x) = E 41 (x-2), (62) 

respectively, o 

We leave open the questions: Are there primes for which there is no 
optimal quadratic polynomial? Are p = 31, 37 such cases? 

4.3 Optimal cubic polynomials 

We now investigate cubic polynomials for the prime numbers po = 2, 3, 5, 7, 11, 
13. 

Theorem 4.6. (i) There is an infinity of optimal cubic polynomials for 
the prime numbers 2, 3. For p = 2 : 

d(x) = x 3 + mx 2 + (pi - 3 - m)x + 2, C x (\) = Pl , (63) 

where m is an arbitrary integer and p 1 a prime number. If |4 + 2m — pi\ is 
prime then Ci(x) is bi-optimal. For p = 3 : 

C 2 (x) = x 3 + [i(pa - 3) - Pl ]x 2 + [2pi - -(pa + 5)]x + 3, 

C 2 (l) = Pl , C 2 (2)=p 2 , (64) 

where p±,p 2 are prime and p 2 odd. If \p 2 — 3p 1 + 3|, \3p 2 — 8pi — 6| are prime 
then C 2 {x) is bi-optimal. 

(ii) There are optimal polynomials forp = 5 : 

C 3 (x) = x 3 - x 2 + 2x + 5, C 4 (x) = 2x 3 + Ax 2 - Ax + 5. (65) 

Po = 7 : 

C 5 (x) = x 3 - x 2 + Qx + 7, (66) 

Po = H: 

C 6 (x) = a; 3 + 5a; 2 + 2a; + 11, C 7 (a;) = a; 3 - 4a; 2 + 5a; + 11, (67) 
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Po = 13 : 

C 8 (x) = x 3 - hx 2 + 8x + 13. (68) 

Proof, (i) Ci(l) = pi is readily verified for p = 2 and 6*2(1) = pi, 6*2(2) = p 2 
for j9 = 3. Note that Ci(-l) = 4 + 2m-p 1 , C 2 (-l) = p 2 - 3pi + 3, C 2 (-2) = 
3p 2 - 8pi - 6. 

(ii) Let C(:r) = ax 3 + bx 2 + cx + p . Then 

Pi = a + b + c + po, p 2 = 8a + 46 + 2c + p , 

p 3 = 27a + 96 + 3c + p . (69) 

(ii) For po = 5, choosing a = l,pi — 7,p2 — 13 yields 6 + c= l, 26 + c = 
and therefore Cz(x) yielding the following additional prime values 

6 3 (3) =29, C 3 (4)=61, (70) 

forming altogether a quintet. 

For a = 2,p = 5, p\ = 7, p 2 = 29 we get c = — b and 6 = 4. The 
polynomial C4 yields the additional prime values 

C 4 (3) = 83, C 4 (4) = 181, (71) 

forming a quintuplet. For p — 7, p\ — 13, P2 = 23 we get b = — 1, c = 6 
and C5 yields the additional prime values 

C 5 (3) = 43, C 5 (4) = 79, C 5 (5) = 137, C 5 (6) = 223, (72) 

a septet and optimal. For p = 11, p 1 = 19, p 2 = 43 we get 6 = 5, c = 2 
and Cq yields the prime values 

C 6 (3) = 89, C 6 (4) = 163, C 6 (5) = 311, C 6 (6) = 419, C 6 (7) = 613, 
C 6 (8) = 859, C 6 (9) = 1163, C 6 (10) = 1531, (73) 

forming an optimal 11-plet. It may be verified that C7 forms an 11-plet also. 
Cs(x) generates the optimal 13-plet 

C 8 (0) = 13, C 8 (l) = 17, C 8 (2) = 17, C 8 (3) = 19, C 8 (4) = 29, 
C 8 (5) = 53, C 8 (6) = 97, C 8 (7) = 167, C 8 (8) = 269, C 8 (9) = 409, 
C 8 (10) = 593, C 8 (ll) = 827, C 8 (12) = 1117. o (74) 

Concluding we ask: Are there optimal cubic polynomials for all prime 
numbers p > 13? 
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